"Noisy" spectra, long correlations and intermittency in wave turbulence 
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We study the k-space fluctuations of the waveaction about its mean spectrum in the turbulence 
of dispersive waves. We use a minimal model based on the Random Phase Approximation (RPA) 
and derive evolution equations for the arbitrary-order one-point moments of the wave intensity 
in the wavenumber space. The first equation in this series is the familiar Kinetic Equation for 
the mean waveaction spectrum, whereas the second and higher equations describe the fluctuations 
about this mean spectrum. The fluctuations exhibit a nontrivial dynamics if some long coordinate- 
space correlations are present in the system, as it is the case in typical numerical and laboratory 
experiments. Without such long-range correlations, the fluctuations are trivially fixed at their 
Gaussian values and cannot evolve even if the wavefield itself is non-Gaussian in the coordinate 
space. Unlike the previous approaches based on smooth initial k-space cumulants, the RPA model 
works even for extreme cases where the k-space fluctuations are absent or very large and intermittent. 
We show that any initial non-Gaussianity at small amplitudes propagates without change toward 
the high amplitudes at each fixed wavenumber. At each fixed amplitude, however, the PDF becomes 
Gaussian at large time. 



I. INTRODUCTION 

The concept of Wave Turbulence (WT), which de- 
scribes an ensemble of weakly interacting dispersive 
waves, significantly enhanced our understanding of the 
spectral energy transfer in complex systems like the 
ocean, the atmosphere, or in plasmas [1—5]. This the- 
ory also became a subject of renewed interest recently, 
(see, e.g. [6-9]). Traditionally, WT theory deals with 
derivation and solutions of the Kinetic Equation (KE) 
for the mean waveaction spectrum (see e.g. [1]). How- 
ever, all experimentally or numerically obtained spectra 
are "noisy" , i.e. exhibit k-space fluctuations which con- 
tain a complimentary to the mean spectra information. 
These k-space fluctuations always develop in numerical 
experiments even though, typically, most numerical ex- 
periments (e.g. [6,7]) start with initial waveaction fields in 
the k-space which have random phases but which have no 
amplitude fluctuations. How fast and why do these am- 
plitude fluctuations get developed? Are they a numerical 
artifact or a real physical phenomenon? Are they Gaus- 
sian or some intermittent bursts of Fourier amplitudes 



can be expected? These questions remain unanswered 
because the waveaction fluctuation have not been stud- 
ied before. Such a study involves description of the higher 
one-point moments of the Fourier amplitudes and it will 
be the main focus of the present work. We will show that 
when these one-point moments are not Gaussian the co- 
ordinate space fields are long-correlated. Such fields are 
very common in WT, and the numerical initial conditions 
discussed above is a typical example. Thus, we will have 
to generalize WT to include such long correlated fields. 



II. RANDOM PHASES VS GAUSSIAN FIELDS. 

The random phase approximation (RPA) has been 
popular in WT because it allows a quick derivation of 
KE [1,3]. We will use RPA in this paper because it pro- 
vides a minimal model for for description of the k-space 
fluctuations of the waveaction about its mean spectrum, 
but we will also discuss relation to the approach of [2] 
which does not assume RPA. 

Let us consider a wavefield a(x, t) in a periodic box 1 



1 Periodic box is an essential intermediate step for formulating RPA and for defining the new correlators Mjf later in this 
paper. This is related to the fact that, strictly speaking, the infinite-space Fourier transform is a distribution, rather than a 
smooth function, for the class of functions corresponding to statistically homogeneous fields. The previous theory considered 
a class of correlators which are box-size independent and which could be formally obtained via a direct manipulation with the 
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of volume V and let the Fourier transform of this field be 
dfc. Later, we take the large box limit in order to consider 
homogeneous wave turbulence. Let us write the complex 
function a k as a k = A k ip k where A k is a real positive 
amplitude and ip k is a phase factor which takes values 
on the unit circle centered at zero in the complex plane. 
By definition, RPA for an ensemble of complex fields a k 
means that 

1. The phase factors ip k are uniformly distributed on 
the unit circle in the complex plane and are statis- 
tically independent of each other 

(V>fciV>fc 2 } = 5\, 
where 8\ is the Kronecker symbol. 

2. The phases are statistically independent from the 
amplitudes A k , 

(ip kl A k2 ) = 

Thus, the averaging over the phase and over the 
amplitude statistics can be performed indepen- 
dently. 

3. the fluctuations of the amplitudes A k must also be 
decorrelated at different fc's. 

(A n k A%) = (Al)(A%) (m,n= 1,2,3, ...). 

Properties 2 and 3 have typically not been men- 
tioned explicitly before. The name RPA itself does 
not refer to the amplitudes but to the phases only. 
However, this important assumption about the am- 
plitude statistics has always been made implicitely 
when using RPA, often without even realizing it. 

To illustrate the relation between the random phases 
and Gaussianity, let us consider the fourth-order moment 
for which RPA gives 

(a kl a k2 a k3 a ki ) = n kl n k2 (6lSl + 5\5l) + QhJl5\5\, (1) 

where 

n k = (4) 



is the waveaction spectrum and 

Q k = (Ai)-2(Al) 

is a cumulant coefficient. The last term in this expres- 
sion appears because the phases drop out for ki = k 2 = 
k3 = k 4 and their statistics poses no restriction on the 
value of this correlator at this point. This cumulant part 
of the correlator can be arbitrary for a general random- 
phased field whereas for Gaussian fields Q k must be zero. 
Such a difference between the Gaussian and the random- 
phased fields occurs only at a vanishingly small set of 
modes with ki = k 2 = k 3 = k 4 and it has been typically 
ignored before because its contribution to KE is negligi- 
ble. Therefore, if the mean waveaction spectrum was the 
only thing we were interested in, we could safely ignore 
contributions from all (one-point) moments 

M { k p) = (\a k \ 2 P) (p= 1,2,3,..). 

However, it is precisely moments M^f* that contain in- 
formation about fluctuations of the waveaction about its 
mean spectrum. For example, the standard deviation of 
the waveaction from its mean is 

e fe = ((K| 4 )-(|a fe | 2 ) 2 )V2 = (Mf -nlfl* (2) 

This quantity can be arbitrary for a general random- 
phased field whereas for a Gaussian wave field the fluctu- 
ation level £fc is fixed, = n k . Note that different values 
of moments Mjf can correspond to hugely different typ- 
ical wave field realizations. In particular, if = n p 
then there is no fluctuations and A k is deterministic, 
£ k = 0. For the opposite extreme of large fluctuations 
we would have 3> n p which means that the typical 
realization is sparse in the k-space and is characterized 
by few intermittent peaks of A k and close to zero values 
in between these peaks. Such an information about the 
spectral fluctuations of the waveaction contained in the 
one-point moments is completely erased from the 
multiple-point moments by the random phases and it is 
precisely why these new objects play a crucial role for 
the description of the fluctuations. 

Will the waveaction fluctuations appear if they were 
absent initially? Will they saturate at the Gaussian level 



infinite-box a k 's. 
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S,k = nk or will they keep growing leading to the k-space 
intermittency? To answer these questions, we will use 
RPA to derive and analyze equations for the moments 
Mjf* for arbitrary orders p and thereby describe the sta- 
tistical evolution of the spectral fluctuations. Note that 
RPA, without a stronger Gaussianity assumption, is to- 
tally sufficient for the WT closure at any order. This 
allows us to study wavefields with moments very 
far from their Gaussian values, which may happen, for 
example, because of the choice of initial conditions or a 
non-Gaussianity of the energy source in the system. 

In [2] non-Gaussian fields of a rather different kind 
were considered. Namely, statistically homogeneous wave 
fields were considered in an infinite space which initially 
have decaying correlations in the coordinate space and, 
therefore, smooth cumulants in the k-space, e.g. 

(a kl ak 2 ak,a ki ) = n kl n k2 {8 k k %l + *M + d^lltl 

where C123 is a smooth function of ki,k2,k3 and <5's 
now mean Dirac deltas. On the other hand, by taking 
the large box limit it is easy to see that our expression (1) 
corresponds to a singular cumulant C123 = Qfci/V 5^5^ . 
It tends to zero when the box volume V tends to infinity 
and yet it gives a finite contribution to the waveaction 
fluctuations in this limit. 2 This singular cumulant cor- 
responds to a small component of the wavefield which is 
long-correlated, - the case not covered by the approach 
of [2] . On the other hand, it would be straightforward to 
go beyond our RPA by adding a cumulant part of the ini- 
tial fields which tends to a smooth function of ki , k2 , k3 
in the infinite box limit (like in [2]). However, such cu- 
mulants would give a box-size dependent contribution to 
the waveaction fluctuations which vanishes in the infinite 
box limit (e.g. it would change £ 2 by Ckkk/V)- Thus, in 
large boxes the waveaction fluctuation for the fields with 
smooth cumulants is fixed at the same value as the for 
Gaussian fields, — nk, and introduction of the singu- 
lar cumulant is essential to remove this restriction on the 
level of fluctuations. On the other hand, the smooth part 
of the cumulant has no bearing on the closure ( as shown 
in [2]) and on the large-box fluctuation and, therefore, 



will be omitted in this manuscript for brevity and clarity 
of the analysis. 

III. TIME-SCALE SEPARATION ANALYSIS 

Consider weakly nonlinear dispersive waves in a peri- 
odic box. Here we consider quadratic nonlincarity and 
the linear dispersion relations Uk which allow three- wave 
interactions. Example of such systems include surface 
capillary waves [4] and internal waves in the ocean [9]. 
In Fourier space, the general form for the Hamiltonian 
systems with quadratic nonlinearity looks as follows, 3 

00 00 
H = ^uj n \c n \ 2 + e ( v inCic m c n S l m+n + c.c.) , 

n— 1 l,m,n— 1 

31~C 

ici = — , ci = aie~ tuJl \ 
oci 

oo 
m,n=l 

+2V l ™a n a m e~ l ^ t 6r +n ) , (3) 

where a n — a(k n ) is the complex wave amplitude in the 
interaction representation, k n = 2nn/L, L is the box 
side length, n = (rai,n 2 ) for 2D, orn= (^1,^2,^3) m 
3D, (similar for h and k m ), uj l mn = u: kl - u km -u km and 
uoi = ujki is the wave linear dispersion relation. Here, 
<~ 1 is an interaction coefficient and e is introduced 
as a formal small nonlinearity parameter. 

In order to filter out fast oscillations at the wave pe- 
riod, let us seek for the solution at time T such that 
2-k/uj <C T <C 1/uj£ 2 . The second condition ensures that 
T is a lot less than the nonlinear evolution time. Now let 
us use a perturbation expansion in small e, 

im\ (°) 1 (!) 1 2 (2) 

Substituting this expansion in (3) we get in the zeroth 
order a\°\T) = a;(0), i.e. the zeroth order term is time 



2 Thus, assuming a finite box is an important intermediate step when introducing the relevant to the fluctuations objects like 

Q fc . 

3 We will follow the RPA approach as presented by Galeev and Sagdeev [3] but deal with a slightly more general case where 
the wave field is not restricted by the condition a(k) = a(—k). We will also use elements of the technique and notations of [2]. 
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independent. This corresponds to the fact that the in- 
teraction representation wave amplitudes are constant in 
the linear approximation. For simplicity, we will write 
a [ (0) = a h understanding that a quantity is taken at 
T = if its time argument is not mentioned explicitly. 
The first order is given by 



OO 

; E 

m,n— 1 



v, "Kn n a ' m Cln ^ m n ^ rn + n 



+2V l ™a m a n hZ8T +n ) , (4) 



where A l mn = e^dt = (e**~» T - l)/iu l mn . Here 



we have taken into account that a|°'(T) = a; and 

4 1} (0) = 0. 

To calculate the second iterate, write 

m,n— 1 



We now have to substitute (4) into (5) and integrate over 
time to obtain 



(5) 



+2V% (V^a m a,aj; + „E[-Lo^,-co™] + 2V£,a m a ll a v E[wg n , <5 ; ™„] , 



m+n 



(6) 



where we used aj^(0) = and introduced E(x,y) — 
J T A(x-y)e i y t dt. 



IV. STATISTICAL DESCRIPTION 

Let us now develop a statistical description applying 
RPA to the fields Since phases and the amplitudes 



are statistically independent in RPA, we will first per- 
form average over the random phases (denoted as (...),/,) 
and then we average over amplitudes (denoted as (...)a) 
to calculate the moments, 

Mjf\T)^(\a k (T)\ 2 ^ A . p=l, 2,3..., 

First, let us calculate \ai(T)\ 2p as 



\ a i( T )\ 2p = (4 0) +«PWaPy (a^+e^+e 2 ^ = \a?^ + e P \a?^-> («<°> a™ + W) 



+ 



e 2 \ ai \ 2 ^ 



where C 2 is the binomial coefficient. 

Up to the second power in e terms, we have 



2p _ 



e 2 \ ai \ 2 P- 2 



(\ai(T)\ 2p U 



(7) 



Here, the terms proportional to e dropped out after the 
phase averaging. Further, we assume that there is no cou- 
pling to the k — mode, i.e. V*j£ = V£* k=0 = 0, so that 

there is no contribution of the term like ((a.^aj 1 ^) 2 )^,. 
We now use (4) and (6) and the averaging over the phases 
to obtain 
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<h (1) l 2 >^ = 

CO 

2£[|v;Ll 2 ^JAU 2 K l | 2 KI 2 

m.n 

+2\Vr m m + S^M 2 \a m \ 2 }, 

oo 

-8|«/| 2 E[l^| 2 ^+^(0,^„)|a m | 2 

m.n 

+l^| 2 <5? +m ^(o, w rj(|aml 2 -KI 2 )]. 

Let us substitute these expressions into (7), perform am- 
plitude averaging, take the large box limit 4 and then large 
T limit (T > l/ujf. We have 



Mi P \T) = M<?\o) +T (-p lk Mi P) +P 2 VkMt 1) ) 



(8) 



with 



Vk = 47re 2 J dkidkaninad^l 2 ^^^) 

+2|^ 1 | 2 ^(u, 2 1 )), (9) 
Ik = Sne 2 J d^dk^V^ 2 5^5^)712 

+ \V k \\Ht 1 S(0(ni-n 2 )). (10) 

Now, assuming that T is a lot less than the nonlinear 
time (T <C 1/W 2 ) we finally arrive at our main result, 



M 



(p) 



-PTfcM^+AfcM^- 1 ). 



(11) 



In particular, for the waveaction spectrum = n k 

(11) gives the familiar kinetic equation (KE) 



n k = -Iktik +Vk = e 2 J(n k ), 
where J(n k ) is the "collision" term [1,3], 

J(n k ) = J dk 2 dk ± (R k 12 - R\ 2 - R\ k ), 



(12) 



with 



Rki2 =^ir\Vf 2 \ 2 5^ 2 8{^ 2 )^n 2 n 1 -n k {n 2 +n 1 )y (13) 
The second equation in the series (11) allows to obtain 



the r.m.s. £? 



M, 



n k (2) of the fluctuations of the 



waveaction (la/^ 2 ). We emphasize that (11) is valid even 
for strongly intermittent fields with big fluctuations. 



V. ANALYSIS OF SOLUTIONS: GAUSSIANITY 
VS INTERMITTENCY 



Let us now consider the stationary solution of (11), 
M^jf 1 = for all p. Then for p = 1 from (12) we have 
Vk = 1 n k- Substituting this into (11) we have 



M 



(p) 



-.pMjf-^nk, 



with the solution M k P ^ = p\ n p k . Such a set of moments 
correspond to a Gaussian wavefield a k . To see how such a 

Gaussian steady state forms in time, let us rewrite (11) in 

(v) 

terms of relative deviations of M k from their Gaussian 
values, 



M 



(p) 



■p\ n\ 



p\ n\ 



p= 1,2,3,.... 



By definition, is always zero. For p = 2, this expres- 
sion measures the flatness of the distribution of Fourier 
amplitudes at each k. This quantity determines the r.m.s. 
^ = — n k of the fluctuations of the waveaction 

(|afc| 2 ), (2) or the mean level of "noisyness", 



£ 2 =n 2 (2^f + 1). 



Using (11), we obtain 



r k 



PVk / p(p-l) 

n k k 



(14) 



4 The large box limit implies that sums will be replaced with integrals, the Kronecker deltas will be replaced with Dirac's deltas, 
5 l m+n — > 5 l mn /V, where we introduced short-hand notation, S l mn — 5(ki — k m — k„). Further we redefine M^f 1 /V p —* Mjf\ 
5 Note that lim E(0,x) = T{tt8{x) + iP(±)), and lim |A(ar)| 2 = 2%TS(x) (see e.g. [2]). 
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for p = 2,3,4, ... This results has a particularly simple 
form of a decoupled equation for p = 2, 

^(2) = 2rjk F (2) 
k n k k 

Taking into account that rj k > 0, we see from this equa- 
tion that deviations of the mean level of fluctuations from 
Gaussianity always decay. In fact, deviations decay 
at each fixed p. This is easy to see from the general solu- 
tion of (14) (obtained recursively): 



-p() 



e p - j p\ 



E 



F { k j) (t = 0), 



(15) 



where 8 — f ^-dt' is a "renormalized" time variable. 

JO n k 

One can see that this expression decays exponentially as 
t — > oo for any fixed p. 

However, an interesting picture emerges at high p cor- 
responding to high wave amplitudes. Although the devi- 
ations F k P ^ eventually decay at each fixed p, their initial 
values propagate in p without decay toward the larger 
values of p. Indeed, one can approximate (14) for p 3> 1 
by a first-order PDE, 



dtF k 



(i>) 



+ ™d p Fj; P) = 0. 
n k p k 



According to this equation, F^ propagates toward high 
p's as a wave. This wave does not change shape with re- 
spect to coordinate x = hip and, therefore, it spreads in 
p without change in amplitude. The speed of this wave 
(in x) is time independent for statistically steady states 
(i.e. when n k and r\ k are time independent). Note that 
this dynamics occurs at each k practically independently, 
i.e. the only coupling of different k's occurs in the prop- 
agation speed via r\ k . 

These solutions allow us to establish the character of 
intcrmittcncy in wave turbulence systems, i.e. to describe 
how high-amplitude "bursts" occur with greater than 
Gaussian probabilities. In terms of the PDF, the wave of 
non-Gaussianity F^ toward high values of p corresponds 
to a wave propagating from low-amplitude "bulk" part 
to the high-amplitude "tail" on the PDF profile. Indeed, 
a Gaussian PDF for a k corresponds to a distribution of 
A = |a fc | 2 of form P(X) = n^e"*/™, and moment M { k p) 
"probes" this distribution in a range of A around A p = pn 
with a characteristic width <5A <~ n. Lifting P(X) in this 



range by a certain factor will result in an increase of mo- 
ment Afj^ by the same factor. Thus, the wave propagat- 
ing from small to large p's corresponds to a wave from low 
to high A's. This wave is such that the relative deviation 
from distribution P(X) = n _1 e _A /™ remains unchanged, 
but the range of A's at which such non-Gaussianity occurs 
moves into the tail (with speed rj) and spreads (propor- 
tionally to its position in A). Note however that at each 
fixed X deviations from P{X) = rT 1 e~ x / n decay, which 
corresponds to decay of F^ at each fixed p at large time. 

Predictions (15) about the behavior of fluctuations of 
the waveaction spectra can be tested by modern experi- 
mental techniques which allow to produce surface water 
waves with random phases and a prescribed shape of the 
amplitude \a k \ [11]. It is even easier to test (15) numeri- 
cally. Consider for example capillary waves on deep wa- 
ter. If a Gaussian forcing at low k values is present, the 
steady state solution of the kinetic equation corresponds 
to the Zakharov-Filoncnko (ZF) spectrum of Kolmogorov 
type [1,4]. It is given by 

n k = Ak- 17 '\ (16) 

with A = \ r Pp 3 / 2 C/<T 1 / 4 , where P is the value of flux of 
energy toward high wavenumbers, p and a are the density 
and surface tension of water, and C ~ 13.98. The sim- 
plest experiment would be to start with a zero-fluctuation 
(deterministic) spectrum and to compare the fluctuation 
growth with the predictions of (11). Note that such no- 
fluctuations initial conditions were used in [6,7]. 

Let us calculate the rate at which a fluctuations grow 
for such an initial conditions. Since n k and r\ k are time in- 
dependent in this case, we have 9 = n k t/n k = j k t. Thus, 
the only quantity we need to calculate is -f k . Let us take 
into account that the spectrum n k is isotropic, that is it 
depends only on the modulus of the vector, not on its 
directions. We then can perform an angular averaging of 
(10) obtaining 

lk = 8e 2 J dK^S kk \ k2 {\V k 2 \ 2 8{u k 2 )n 2 



S k u = ^(k-ki-ka)) 



+ \V k 2 1 \ 2 S(u,l 1 )(n 1 -n 2 )), 
J 5(k - kx - k 2 ) dBxdBi 



- {{kk^ 2 + (kk 2 ) 2 + (fcife) 2 ) - fc 4 - k\ - k\ . 



(17) 



Let us substitute ZF spectrum (16) into (17), take the 
values of L0 k and V k 2 appropriate for the capillary waves 
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on deep water( [1], eqs (5.2.1-2)). By changing the vari- 
ables of integrations via k\ — fc£i, ki — kt^ we can 
factor out the k dependence of 7^ . Performing one of £ 
integrals analytically with the use of the delta function in 
o/s, we perform the remaining single integral numerically 
to obtain (all the integrals converge): 

4.3(hVa 3/4 
7 16np 3 / 2 ' 

where the dimensionless constant 4.30 was obtained by 
numerical integration. Substituting the value of A we 
finally obtain 

7 = LWVFaVW, 

Consequently, our prediction for the fluctuations growth 
is 

= nl(2F^ + 1) = A 2 k- l7 ' 2 (l - e- 2 ^). (18) 

Note that fluctuations stabilize at Gaussian values 
faster for high k values. One can also substitute 7fe cal- 
culated for the capillary waves into the solutions for the 
higher p's, (15). Again, the dynamics here is going to 
be faster at large fc's because they correspond to higher 
values of 6 = jt. In particular, at large k there will 
be a faster wave toward higher p's. For the particular 
type of initial conditions we have taken (no initial fluctu- 
ations), this wave will describe formation of a Rayleigh 
distribution (corresponding to the Gaussian statistics of 
dfc) behind a propagating front on the PDF profile. It a 
way, this dynamics is non-intermittent: zero initial fluc- 
tuations grow to the Gaussian level but never exceed it. 

It is also interesting to test our predictions when the 
initial conditions, or forcing, are non-Gaussian, as in 
most practical situations. Our theory predicts that non- 
Gaussianity of the low-amplitude (bulk) part of the PDF 
will propagate without decay into the high-amplitude tail 
at each fixed k. The speed of this propagation is propor- 
tional to 7^ and, therefore, will be higher for large fc's in 
the case of the capillary waves. This means higher in- 
termittency in the low- A; range in the case of stationary 
forced turbulence. 

VI. DISCUSSION 

In this manuscript, we derived a hierarchy of equa- 
tions (11) for the one-point moments of the wave- 



action \dk\ 2 ■ This system of equations has a "triangu- 
lar" structure: the time derivative of the p-th moment 
depends only on the moments of order p,p — 1 and 1 
(spectrum). Their evolution is not "slaved" to the spec- 
trum or any other low moments and it depends on the 
initial conditions. RPA allows the initial conditions to 
be far from Gaussian and deviation of n'th moment from 
its Gaussian. Among two allowed extreme limits are the 
wavefield with a deterministic amplitude \a k \ (for which 
M^jf 1 = n p k ) and the intermittent wavefields character- 
ized by sparse k-space distributions of | (for which 
M^»nl). 

Equations (14) for the deviations from Gaussianity 
have an interesting property that the nonlinear coupling 
between different modes k occurs only via a rate constant 
rj/n. By removing this dependence into a "renormahzed" 
time 9 one gets a linear system of equation which can be 
easily solved in the general case, see (15). Analyzing 
these solutions we showed that the deviation from Gaus- 
sianity decreases as at each fixed amplitude \ak\- At the 
same time, we showed that any initial non-Gaussianity 
at small amplitudes propagates as a non-decaying wave 
toward the high-amplitude tail of the PDF. This pro- 
cess describes the character of the wave turbulence inter- 
mittency when high-amplitude wave "bursts" occur in 
the system more frequently than predicted for Gaussian 
fields. On the other hand, the assumption about the weak 
nonlincarity breaks down when a high amplitude burst 
occurs in the system, leading to a failure of the RPA clo- 
sure to describe the PDF tails. One can conjecture that 
the resulting phase coherence will lead to a nonlinear am- 
plitude saturation which will stop the wave predicted by 
our theory which, in turn, will lead to a stagnation and 
accumulation in this region on the PDF tail. Thus, it is 
natural to expect even stronger intcrmittency when the 
higher order nonlinear effects are taken into account. 

We would like to emphasize that the type of intermit- 
tency discussed in the present manuscript appears within 
the weakly nonlinear closure and not as a result of its 
breakdown as in [10]. This intermittency is quite subtle 
and it occurs only in the PDF tails and not in its core 
(which tends to a Gaussian state). As a result, the lower 
moments will not feel these rare "bursts" and they will 
evolve as predicted by the WT closure. We have showed 
that this kind of intermittency inevitably occurs at all 
wavenumbers k provided some initial non-Gaussianity is 
present in the PDF core. Paper [10] considers a different 
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and a more dramatic kind of intermittency which occurs 
simultaneously with the strong nonlinearity of the typi- 
cal wave from the PDF core. This kind of intermittency 
is more seldom and it takes place only in some special 
parts of the fc-space (e.g. at very small scales). In par- 
ticular, it never occurs for the capillary waves considered 
in this paper provided that only weakly nonlinear waves 
are produced at the forcing scale. 

The present paper deals with the three-wave systems 
only. The four-wave resonant interactions are slightly 
more complicated in that the nonlinear frequency shift 
occurs at a lower order in nonlinearity parameter than 
the nonlinear evolution of the wave amplitudes. To build 
a consistent description of the amplitude moments one 
has to perform a renormalization of the perturbation se- 
ries taking into account the nonlinear frequency shift. 
This will be done in a future publication. 
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